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VALIDATION OF A PSEUDO-SOUND THEORY FOR THE PRESSURE-DILATATION 
IN DNS OF COMPRESSIBLE TURBULENCE 


J. R. RISTORCELLI* AND G. A. BLAISDELL^ 

Abstract. The results of an asymptotic theory for statistical closures for compressible turbulence are 
explored and validated with the direct numerical simulation of the isotropic decay and the homogeneous shear. 
An excellent collapse of the data is seen. The slow portion is found to scale, as predicted by the theory, with 
the quantity Af 2 and e 9 . The rapid portion has an unambiguous scaling with a 2 M 2 e s [^- — 1](§^) 2 - Implicit 
in the scaling is a dependence, as has been noted by others, on the gradient Mach number. A new feature 
of the effects of compressibility, that of the Kolmogorov scaling coefficient, a, is discussed. It is suggested 
that a may contain flow specific physics associated with large scales that might provide further insight into 
the structural effects of compressibility. 
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1. Introduction. The validation, using DNS data, of the scaling predicted by an analytical develop- 
ment for the pressure dilatation, < pd >, appearing in the single point closures for compressible turbulence, 
given in Ristorcelli (1995, 1997), is the subject of this article. The analysis is relevant to shear flows with 
negligible bulk dilatation and low M t 2 . These restrictions are satisfied in a wide number of flows ranging 
from simple shear layers of theoretical interest, Papamoschou and Roshko (1988), to the complex shear layers 
associated with supersonic mixing enhancement, Gutmark et al. (1995). In most of these supersonic shear 
layers a Mach number based on the fluctuating velocity of the fluid particle is small. A Mach 4 mean flow 
with a turbulence intensity of 8 per cent has a turbulent Mach number of M t = 0.32. The square of this 
turbulent Mach number, the appropriate perturbation expansion parameter arising from the Navier-Stokes 
equations, M 2 ~ 0.1, is small. The existence of this small parameter, Af 2 , allows some analytical results. 

The article is primarily a study of the analysis of Ristorcelli (1995, 1997) in the light of recent DNS 
results of Blaisdell. The representations given in Ristorcelli (1995, 1997) were obtained using scaling argu- 
ments about the effects of compressibility, a singular perturbation idea and the methods of statistical fluid 
mechanics. While the results are expressed in the context of a statistical turbulence closure, they provide, 
with few phenomenological assumptions, an interesting and clear physical model for the scalar effects of com- 
pressibility. For a homogeneous turbulence with quasi-normal isotropic large scales the expressions derived 
are - in the small turbulent Mach number squared limit - exact. The analytical results, which are a rigorous 
consequence of the low Af 2 assumptions and do not contain any unspecified empirical coefficients, are shown 
to predict the scalings in the DNS of homogeneous compressible turbulence. 
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2. Preliminaries. It has been shown, Ristorcelli (1995, 1997), that a low fluctuating Mach number 
expansion for the compressible Navier Stokes equations produces a diagnostic relationship for the dilatation. 
The small parameter in these expansions is the square of the fluctuating turbulent Mach number: 7 Mf and 
M t — u/coo, where u = 2k/3 = {ujUj )/ 3 and = 7 Poo/ Poo- To leading order, the density fluctuations are 
given, in nondimensional units, by the solenoidal pressure fluctuations, 7/01 = p inc = p x and the continuity 
equation becomes a diagnostic relation for the fluctuating dilatation, 

(!) -id - p, t +v k p,h • 


The subscript on p\ has been dropped. It is seen that one does not need to obtain a solution to the evolution 
equation for the compressible velocity field, uij, in order to obtain its dilatation, d = A very useful result. 
The dilatation is diagnostically related to the local fluctuations of the pressure and velocity; it is the rate of 
change of the incompressible pressure field p U jj = (viVj)^, following a fluid particle. Constitutive relations 
for the pressure-dilatation can be found by taking the appropriate moments of the diagnostic relation for 
the dilatation to produce, 

(2) -27 < pd > =^ t <pp> . 

The near field compressibility effects, as manifested in < pd >, have been directly linked to the solenoidal 
portions of the velocity field. This fact has been exploited to obtain expressions for the pressure- dilatation 
covariances, Ristorcelli (1995, 1997): 


( 3 ) 


where 


( 4 ) 


<pd>= - XpdM? [P fe -pe + T k ] - pk M 2 Xp d ~T 
+ ^X P d^t 4 7(7 - l)(Pr + pe + T t )] 

2Ipd 

Xpd 1 + 2 I^M? + § 7 ^ 7(7 - 1 ) 

x r d = ^ 

P 1 + 2 I pd M 2 + § IpdM^ - 1) 

Ipd = + 7 ^ [3S 2 + 5W 2 ] 

jr __ ]L(-\3 2 jr 

pd 30 3 7l * 


Note that T = [3*S 2 + 5W 2 }. The nondimensional strain and rotation rates are given by: S 2 = ( Sk/e s ) 2 , 
IV 2 = (Wk/e s ) 2 where S = yj S lJ S lJ and W = y/W tJ W tJ . The strain and rotation tensors are defined 
in analogy with the incompressible case, i.e., traceless S {j = W i3 = -U hi \. 

Here M t — u/c where u = 2k/3 =< UjUj > /3, = l Poo/ Poo and a comes from the Kolmogorov 

scaling relation £ = a(2fc/3) 3 / 2 /e 9 . The constants, denoted by the /*, in these expressions are given by 
integrals of the longitudinal correlation function, /: If = / 0 °° and If = 2 / 0 °°f/d£. A quick order 

of magnitude estimate for the integrals can be made using / — *7 4 The following values are found: 

If — ~ = 1.273. The values found from high Reynolds number wind tunnel data are similar: 
If = 0.300, If = 1.392, Zhou (1995). 
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3. Isotropic decay. For the isotropic decay the expression for the pressure-dilatation becomes 

( 5 ) <pd>= XpdM f e 8 . 

Here p has been set to unity. The sign of < pd > is positive indicating a net tranfer of energy from potential 
to kinetic modes. Rearranging, to isolate the scaling, produces 

(6) < pd > = x , = ill 

M?e s Xpd 1 + | I* Ml 

Terms of order Mf have been dropped. Earher estimates given in Ristorcelli (1995), shown above, indicate 
If = 0.5 - 0.3. The theory therefore predicts an asymptote for Xpd as the turbulent Mach number vanishes: 

(7) Xpd -> 0.666 - 0.40 as M 2 -> 0. 

The DNS results, shown in Figure 1, were provided by Blaisdell for three different initial turbulent Mach 
numbers. As a service to the reader, the figure identifies two definitions of the turbulent Mach number: that 
used by Blaisdell et ai (1993) in his simulations, and that which comes from the perturbation analysis 
of Ristorcelli (1995, 1997). The present compressible DNS reflect a consistent set of initial conditions as 
described in Ristorcelli and Blaisdell (1997). 

The agreement with the DNS, shown in Figure 1, is very good. The theory has been corroborated without 
a posteriori adjustment of constants. The actual values of the constant could in principle be calculated from 
the DNS. As they are expected to be weakly dependent on initial conditions, this course is not followed 
further - what has been presented is sufficient for verification. Moreover, the slow portion of the pressure- 
dilatation is nominal compared to the rapid portion which is the most important contributor in the shear 
flows of interest. 


4. Homogeneous shear. The pressure-dilatation in the homogeneous shear is now investigated. The 
instantaneous pressure- dilatation is seen in Figure 2. Also shown are its averaged values following the 
procedure of Sarkar (1992). Here the time integral of the pressure-dilatation has been taken: the vertical 
axis being J ^^d(St). The oscillations in the pressure-dilatation associated with the relaxation from 
initial conditions are not seen. There is, nonetheless, a build up of the oscillations which has been linked to 
the compressible component of the pressure field, Sarkar (1992), Blaisdell and Sarkar (1993). An explanation 
consistent with this observed behavior has been advanced in Ristorcelli (1997b). 

For homogeneous shear, the expression for the pressure-dilatation can be simplified. For simple shear 
T = 8S 2 and, neglecting terms of order M 4 , one obtains 

(8) <pd>= —\pdM? [P k -e] -k M? X r pd 8§~ t S 2 . 

Here S — For Blaisdell’s homogeneous shear, in which S = const . , the expression can be rearranged 


(9) 


<pd >= - XpdM f e 9 [— - 1] 

- 16 xU M t £ s (— ) 2 D 
£ 


Sk 


D(St) £ s * 

Note that the coefficent of the first term, Xpdi ignoring the small slow pressure contribution, scales as 
Xpd ~ (ff) 2 ; accounting for the definitions of the x’s, the pressure-dilatation can be rewritten as 


( 10 ) 


< pd - 1)1 ^ P+ a 1 _ lC( M) fl 
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As the flow evolves, it is expected that I{ — ► const and [1 + -p~ — r>(st) f^] 1* The scaling of < pd > with 

the term in the first set of brackets will be investigated. 

4.1. “Non-equilibrium” aspects of < pd >. The homogeneous shear DNS is not an equilibrium 
flow. Using as an indicator of the non- equilibrium nature of the flow one sees that 7^ — > 0 only in 
the latter stages of the DNS, see Figure 3 bottom plot. Two regions of the flows evolution are accordingly 
distinguished: a nonequilibrium earlier portion in which D/Dt S ^ 0 , and a structural equilibrium portion 
for which D/Dt S — ► 0. If the flow is in structural equilibrium, ^ — ► const , one obtains the simplest form 
of the pressure-dilatation model: denote it by the subscript “se” , 

(11) < pd >se= ~ XpdM 2 (Pfc - e] . 


Let this be called the structural equilibrium form of the pressure-dilatation (which does not mean that the 
flow is in equilibrium as Pk ^ s). 

Figure 3 indicates the relative contributions of the two terms making up the the pressure- dilatation 
model. It is seen at small times, that both terms make non-negligible contributions to the pressure- dilatation. 
As the structural equilibrium is approached, D/Dt S —+ 0 , and the second term’s importance, as might be 
expected, becomes negligible. This is manifested in the second graph were the ratio of < pd > se to the total 
< pd > approaches unity. Also shown is the percentage time rate of change of the relative strain §£ 757 • 

The curves are noisy as they involve differentiation of numerical data; the trends are nonetheless apparent. 


4.2. Pressure-dilatation scalings. The appropriately scaled integrals of < pd > will now be taken. 
In this way one can establish whether the scalings predicted by the model are correct. In the latter portions 
of Blaisdell’s DNS, St > 10 , about three to four eddy turnovers past initialization, a structural equilibrium 
is approached. In this region the scaling suggested by 


( 12 ) 

is investigated. 
J ST denote this 


<pd> [a 2 ( — ) 2 M 2 £,(— - 1 )] /[ 

After the scaling < pd > the time integral /( )d(St) will be taken. Let the symbol 

averaging operation. The integrals 


I 0 


h 


h 


h 


h 


Is 

Is 

Is 

Is 

Is 


ST 


ST 


< pd> 
<pd> 

£3 

<pd> 


[^- 1 ] 

< pd> 


ST M 2 £ a [^-l](f)2 
< pd > 

ST Q 2A/2 £ ,[T t -l](f)2 


will be computed from the data. The integrals are shown in Figure 4. The integration starts at St = 9; the 
integrals are normalized by their values at St = 10 . If the scaling suggested by the analysis of Ristorcelli 
(1995, 1997) is correct the last integral, J 4 , will be approximately constant, reflecting the fact that the time 
integral ^ f 7[d(5t) ~ I{ and -^S — » 0 as the equilibrium portions of the DNS are attained. Inspection of 
I\ in Figure 4 shows this to be the case. The period of time 10 < St < 16 corresponds to about one eddy 
turnover time, k/e. 
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4.3. The gradient Mach number. The largest relative collapse of the scaled integrals of < pd > 

occurs with the quantities M t 2 (^) 2 - the collapse from J 2 to J 3 . The quantity M t 2 (^) 2 can be interpreted 

of as the square of the gradient Mach number. The pressure-dilatation is a strong function of the gradient 

Mach number. Sarkar (1995) has defined a gradient Mach number as M g = S£/c ; the transverse two-point 

correlation of the longitudinal velocity is used as the length scale, t. In this article £ will be taken as the 

traditionally defined longitudinal length scale that occurs in the Kolmogorov scaling: l = a(2fc/3) 3 / 2 /e s . In 

which case a mean strain Mach number is defined: 

/1oX S£ 2Sk 2- 2 - 

(13) Ms — — = a- — M t = a-SM t £* -SM t . 

C O £ ^ O O 

In fact, the curve overshoots the optimum collapse (a horizontal line): the gradient Mach number is increasing 
faster than < pd >. 

4.4. The Kolmogorov scaling coefficient. The collapse, / 4 , is much better when the Kolmogorov 
coefficient is included. A new feature associated with compressibility, that of the Kolmogorov scaling co- 
efficient, a, is thus apparent. The values used for a come from Blaisdell’s DNS. The longitudinal integral 
length scale, £, k and e s are calculated from the DNS. The Kolmogorov relationship, e — a(| k)*/ 2 /£ is used 
to find a. 

The pressure- dilatation model is sensitive to variations in the Kolmogorov scaling parameter: in fact 
it is decreasing as rapidly as Ms is increasing. The definition for the mean gradient Mach number given 
above, Ms = |~M*, implicitly assumes a ~ 1. M t and ~ are not new quantities for describing turbulence 
in single point closures; a, however, is new. This distinction is made in order to recognize a as a new 
independent quantity. 

The Kolmogorov constant is thought to be a universal constant for high Reynolds number in isotropic 
turbulence; for non-ideal, finite Reynolds number, anisotropic turbulence a is a flow specific quantity. It is 
this fact that makes developing a turbulence model from the analysis of Ristorcelli (1995, 1997) difficult: a 
choice for a must be made and for any given flow the choice is not, a priori , known. The value of a can be 
thought of as describing some large scale structural aspects of the flow: it, after all, relates the kinetic energy, 
its cascade rate and the two-point correlations. These ideas are further developed in Ristorcelli (1997c). 

5. Summary and Conclusions. The pressure-dilatation is found to be a nonequilibrium phenomena. 
It scales as a 2 M?(^) 2 {P k /e s — 1]. For it to be important requires both 1) the square of the gradient Mach 
number, A/J, to be substantial and 2) for the flow to be out of equilibrium P k ^ e. The pressure- dilatation 
has been observed to be either positive or negative. Its dependence on the non-equilibrium nature of the flow, 
Pk 7^ £> indicates that the pressure dilatation can be either a stabilizing or destabilizing. These predictions 
are consistent with the DNS of Simone et ai (1997) who observes such behavior as related to the anisotropy, 
fri2 = 

Except for the well-established Kolmogorov scaling, £ ~ (|fc) 3 / 2 /£ a , and the quasi-normal assumption, no 
additional phenomenological assumptions are made. The pressure- dilatation is a function of the Kolmogorov 
scaling coefficient and this is expected to be an important feature in models for the effects of compressibility 
on turbulence. The Kolmogorov coefficient is a flow dependent quantity: there is little known about its 
dependence in non-ideal anisotropic, strained, inhomogeneous — flow situations. The appearance of the 
Kolmogorov coefficient, in as much as it links the energy, the spectral flux and a two-point length scale, is 
an indication of dependence on large scale structure. 

In all likelyhood, there are many aspects of compressibility that will contribute to the unusual and 
dramatic stabilizing effects of compressibility. It is unlikely that any one term in a statistical closure will 
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account fully for the diverse physical mechanisms. This article has focused on one of these effects, the 
pressure-dilatation. The present analysis treats only the “scalar” effects of compressibility - the reduction k 
through the dilatational covariances in the energy budget; it cannot account for the reduction in the shear 
anisotropy, 612, or the normal anisotropy, 622, so important to the production mechanism for the shear stress, 
< V\V 2 >• To account for these more substantial structural effects appears to require a compressible pressure- 
strain representation accounting for the effects of compressibility. This has been indicated in Blaisdell and 
Sarkar (1992), Vreman et al. (1996) and also Simone et al. (1997). 
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Fig. 3. Contribution of the equilibrium portion of the pressure- dilatation to the pressure- dilatation as predicted by the 
pseudo-sound analysis. 
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